High-momentum dynamic structure function of liquid 3 He— 4 He mixtures: a 
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The high-momentum dynamic structure function of liquid 3 He- 4 He mixtures has been studied 
introducing final state effects. Corrections to the impulse approximation have been included using 



a generalized Gersch- Rodriguez theory that properly takes into account the Fermi statistics of 3 He 



o 

^ . atoms. The microscopic inputs, as the momentum distributions and the two-body density matrices, 

correspond to a variational (fermi)-hypernetted chain calculation. The agreement with experimental 
data obtained at q = 23.1 A -1 is not completely satisfactory, the comparison being difficult due 
to inconsistencies present in the scattering measurements. The significant differences between the 
experimental determinations of the 4 He condensate fraction and the 3 He kinetic energy, and the 
& , theoretical results, still remain unsolved. 

a 

g 

o 

o 



PACS numbers:67.60.-g, 61.12.Bt 



I. INTRODUCTION 



Liquid 3 r He.- 4 Hc mixtures at low temperature have been of long-standing interest from both experimentafli and 
theoreticalErO viewpoints. From the theoretical side, isotopic 3 He- 4 He mixtures manifest fascinating properties intrin- 
sically related to their quantum nature. The different quantum statistics of 4 He (boson) and 3 He (fermion) appear 
reflected in the macroscopic properties of the mixture as its very own existence in the zero-temperature limit. One 
^sO I of the most relevant features that the 3 He- 4 He mixture shows is the interplay between both statistics driven by the 
C*~) 1 correlations. Signatures of that are, on the one hand, the influence of 3 He on the condensate fraction (no) and the 
superfluid fraction (p s / p) of 4 He, and on the other, the change in the momentum distribution of 3 He atoms due to 
the correlations with 4 He. Both theory and experiment show that the 4 He superfluid fraction decreases with the 3 He 
concentration [x) in .the mixture&Q whereas the condensate fraction no moves on the opposite direction showing an 
enhancement with xB Concerning the 3 He momentum distribution in the mixture, microscopic calculationsQ't3 point 
to a sizeable decrease in the values of n(k = 0) and Z = n(kp) — n(kp) with respect to pure 3 He, with a subsequent 

S population at high k. The long tail of n(fc) gives rise to a 3 He kinetic energy which is appreciably larger than in pure 
Mle. 

Experimental information on the momentum distributiorL_tt(£;) can be drawn from deep inelastic neutron scattering 
(DINS), tj £3 as first proposed by Hohenberg and Platzman.tJ It is nowadays well established that at high momentum 
transfer q the scattering is completely incoherent and accurately described by the impulse approximation (IA). As- 
suming IA, the momentum distribution can be directly extracted from experimental data. However, this procedure is 
not straightforward because of the unavoidable instrumental resolution effects (IRE) and the non-negligible final-state 
effects (FSE). The FSE are corrections to IA that take into account correlations between the struck atom and the 
medium which are completely neglected in the IA. At the typical values of the momentum transfer used in DINS on 
helium (q ~ 20 A -1 ), both IRE and FSE broaden significantly the IA prediction hindering a neat detpniiaation of 
n(k). The dominant contributions to the FSE are well accounted by the different theoretical methodscJlij used in 
their study with an overall agreement for q > 16 A _1 Ji3~li2 Using the theoretical prediction for the FSE and the IRE 
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associated to the precision of the measurements, DINS in superfiuid 4 He points to a condensate fraction no = 9.2 ±1.1 
and a single-particle kinetic energy T/N = 14.5 ± 0.5 K.EHl Beit values are in a nice agreementp-with the theoretical 
values obtained withjGreen's function Monte Carlo (GFMC)J^I diffusion Monte Carlo (DMC)E3 and path integral 
Monte Carlo (PIMC)B methods. „ 

Normal liquid 3 He has also been studied by DINS.cil This system is more involved from a technical point of view 
due to the large neutron absorption cross section of 3 He atoms which significantly reduces the signal-to-noisp^atio 
of the data. Recent measurements of S(q,u>) at high q point to a single-particle kinetic energy of 10 ± 2 Kp3'E-3 a 
value that is clearly larger than a previous DINS determination (8.1 ± 1.7 K)Bj A recent theoretical determination 
of T/N using DMC predictsc3 a value 12.24 ± 0.03 K in close agreement with other microscopic calculations 
Therefore, theory and experiment have become closer but the agreement is still not so satisfactory as in liquid 4 He. 
These discrepancies have been generally attributed to high-energy tails in S(q, u>), that are masked rby the background 
noise, or even to inadequacies of the gaussian models used to extract the momentum distributionr^.. . 

In recent years, there have been a few experimental studies of liquid 3 Hc- 4 He mixtures using DINS.EE3 The response 
of the mixture has been measured at two momenta, q = 23.1 A" 1 and q — 110 A -1 , and different 3 He concentrations. 
By using the methodology employed in the analysis of the response of pure phases, results for the He condensate 
fraction and the kinetic energies of the two species are extracted as a function of x. This analysis points to a 
surprising result of no = 0.18,3 a factor two larger than in pure 4 He. Concerning the kinetic energies, a remarkable 
difference between 4 He and 3 He appears. The 4 He kinetic energy decreaseSplinearly with x whereas 3 He atoms 
show a x-independent kinetic energy that is the same of the pure 3 He phaseE'EHl Except for the 4 He kinetic energy, 
those experimental measurements yield values that are sizeably different from the available theoretical calculations. 
Microscopic approaches to the mixture using both variational hypernetted-chain theoryfl (HNC) and diffusion Monte 
Carldij point to a much smaller enhancement of n in the mixture, and a 3 He kinetic energy much larger at small x 
and that decrease with x down to the pure 3 He result. 

The theoretical estimation of the FSE in the scattering is of fundamental interest, and also unavoidable u3-,the 
analysis of the experimental data. In a previous work,tZl we recovered the Gersch-Rodriguez (GR) formalismtfJ for 
liquid 4 He, and proved that using accurate approximations forp-the two-body density matrix the FSE correcting 
function is very close to the predictions of other FSE theories .Il3ll3 The generalization of this theory to a Fermi 
system as liquid 3 He is not straightforward. The convolutive scheme developed for a boson fluid is now impeded by 
the zeros present in the Fermi one-body density matrix. Recently, it has been proposed an approximate GR-FSE 
theory that incorporates the leading exchange contributions without the aforementioned problems.^ 1 ] That theory is 
expected to capture the essential contributions of Fermi statistics, and thereby to be accurate enough to generate the 
FSE correcting function in dilute 3 He- 4 He liquid mixtures. The aim of the present work is to provide microscopic 
results on the FSE effects in 3 He- 4 He mixtures. The inclusion of FSE on top of the impulse approximation allows 
for a reliable prediction on the dynamic structure function at high momentum transfer that can be compared with 
scattering data. 

In the next section, the FSE formalism for the mixture is presented. Sect. Ill is devoted to the lowest energy- 
weighted sum rules of the response and the FSE correcting functions. The results and a comparison with available 
experimental data are reported in Sect. IV. A brief summary and the main conclusions comprise Sect. V. 

II. FSE IN FERMION-BOSON 3 HE— 4 HE LIQUID MIXTURES 

The dynamic structure function of the mixture S(q, lo) is completely incoherent if the momentum transfer is high 
enough. The incoherent total response can be split up in terms of the partial contributions S^ a \q,uj) 

S(q, uj) = 0-4(1 — x) (q, lo) + a 3 x S^{q,w), (1) 

x = N3/ N being the 3 He concentration, and 04, 173 the cross sections of the individual scattering processes (173 = 5.61 
barn, 04 = 1.34 barn). Notice that in this regime the cross term S^ 3 ' 4 \q, to) does not appear because it is fully coherent, 
and the incoherent density- and spin-dependent Fermi responses are identical ( 03 = 03. ^ + a^j with as,d = 4.42 barn 
and &3j = 1.19 barn). In Eq. (|l|), each single term can be obtained as the Fourier transform of the corresponding 
density-density correlation factor 

S( a \q,t) = — J2(e- i,l ^e mt e ici ^e- im ) , (2) 
Na 3=1 

where a =3, 4 stands for 3 He and 4 He, respectively. In Eq. (^), H is the hamiltonian of the system (K = 1), 
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N 3 ,Ni 



H 



2ni4 



Ev?-^Ev? + J E E v {a ' 0) (nj) , 



3 = 1 



(3) 



a,/3=3,4 i,j'=l 



with V^'^inj) the pair-wise interatomic potentials, that in an isotopic mixture, as the present one, are all identical. 
The translation operators act on the hamiltonian H and transforms Eq. (Q) into 



(4) 



with Wq = q 2 1 (2m a ), and Ly' — v*-"' • pj being the projection of the momentum of particle j along the direction 
of the recoiling velocity v( Q ) = q/rn a . One can then define a new operator C^ a \t) 



CW(i) = e -iHt e i(H + L^)t e -iLi"h (5) 

that contains the FSE corrections to the IA. In terms of C a (t), the density-density correlation factor turns out to be 

S (Q %t) = e^^C^itle" 1 "') . (6) 

In the high-momentum transfer regime, in which we are interested in, v^ Q ' is large while t is short, in such a way 
that their product s = v( Q 'i is of order one. In terms of this new variable s, Eqs. (pi) and df!) become 



(7) 



C {a) {s) = e 



sH is(«+v (o, -pi) -isv^-pi 



with a new hamiltonian H = H/v^ a \ The operators C^ a \s) satisfy the differential equation 

d 



d,s 



with 



C (a) (.s) = i[A t (s)WA(a) - H] C {a \s) 



At(s) = e --w e -(«+v ( "»-P 1 ) 



(8) 



(9) 



(10) 



The differential equation (jTJ) may be solved by means of a cumulant expansion in powers of l/v^^Ell In the high-g 
limit, only the first terms of the resulting series are expected to significantly contribute. In fact, the IA is recovered 
when only the zero-order term is retained, 



(11) 



p < f\s) being the one-body density matrix. By including the next-to-leading term, the leading corrections (FSE) to 
the IA are taken into account, 



(12) 



with Hq(s) = e ls ' Pl He * s ' pi . A Gersch-Rodriguez cumulant expansion of Eq. (12) for the 4 He component in the 
mixture leads to a FSE convolutive scheme 



SW(q,s) = S$(q,s)RW(q,s) 



(13) 



with the IA response ([Tl]), and 
i? (4) (<?,s) =cxp 



p! 4) ( S ) 



dr pj 4 ' 4 ' ( r i 0; r + s) 



I -'^s.l-py / d S 'AV(r,s') 



I— J drp^ 3) (r,0;v + s) 1 - exp f-jL J ds'AV(r,s' 



(14) 
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In the above equation, AV r (r,y,r / ) = V{vij + r') — V{rij). Apart from pf\ R^(q,s) is a function of the (4,4) and 
(4,3) components of the semi-diagonal two-body density matrix 



4 a,/3) (ri,r 2 ;ri,r 2 ) = N a (N - S af3 ) 



J dr N - 2 Wg(n, r 2 , ■ ■ ■ , rjv)^^, r 2 , ...,r N ) 
Jdv* | *(n,r 2 ,...,r^) | 2 



(15) 



The analysis of the 3 He (fermion) component is much more involved. A fully convolutive formalism is now forbidden 
because the zero-order cumulant, which is proportional to the one-body density matrix, has an infinite number of 
nodes. Nevertheless, it is plausible to assume that at high q the FSE are dominated by dynamical correlations, and 
that statistical corrections to a purely FSE scheme can therefore be introduced perturbatively. With this hypothesis, 
the 3 He response can be split up in two terms,! 3 -!! 



using the following identity for the n-body density matrix of the mixture 



p N (r 1 ,r 2 , . . .,r N ;r[) = pf\r lv ) 



1 



PiM 



plr(ri,r 2 , . . . ,r;v;ri) 



(16) 



(17) 



p N (r 1 ,r 2 , . . .,rjv;ri) 



Pi(rii') 



p%(r 1 ,r 2 , ■ ■ .,rjv;ri) 



The superscript B stands for a boson approximation, i.e., a fictitious boson-boson 3 He- 4 He mixture. In that factor- 
ization (p7|), the first term allows for a description of the 3 He response in which the IA is the exact one while the FSE 
are introduced in a boson-boson approximation. Statistical corrections to the FSE are all contained in the second 
term. 

In Eq. (16), S*g (q, s) is the main part of the response and can be written as a convolution product 



Sg\q, S ) = S^(q, s)R {3) (q,s) 



with {Qi s ) 



A 3 h„ „\ _ e( 3 ) 

; ISU « / u< } P^i\s)/p3 the impulse approximation, and 
1 



(18) 



R {3) (q,s) 



exp 



drp 2 3 ' 3)B (r,0;r + s) 



1 — exp 



i 



ds' AV(r,s'] 



1 



dr p 2 3,4 ^ B (r, 0; r + s) 



l-exp( I ds' AV(r,s') 



(19) 



the boson-like FSE correcting function. 

The additive correction AS^ 3 \q, s) in Eq. ([l6|) takes into account the statistical exchange contributions in the FSE 
and is expected to be small. Actually, it is a function of 



P2 



(ri,r 2 ;ri) , 



(20) 



according to the decomposition (Jl 7|) - The variational framework of the (fermi)-hypernetted chain equations (F)HNC 
that is used in this work to calculate the one- and two-body density matrices, provides a diagrammatic expansion 
to estimate Ap 2 3 ' a \ Following the diagrammatic rules of the FHNC/HNC formalism, Ap 2 S '°^ may be written as the 
sum of two terms: 



Ap 2 3 ' Q) (r l5 r 2 ; r[) = p a p^ (r n ,)G^ (r u r 2 ; r[) - PaPmMF^n, r 2 ; t{) 



(21) 



p^\r) is the one-body density matrix and Pid(^) is an auxiliary function, which factorizes in p^\r), and that sums 
up all the diagrams contributing to p{\r) except those where the external points 1 and 1' are statistically linked.0 
F( 3 < Q ) and G< 3 ' Q > in Eq. (§|) sum up diagrams with the external vertices (1,1', 2) with and without statistical lines 
attached to 1 and 1', respectively. With this prescription for Ap 2 3 ' a \ the additive term AS^(q, s) becomes finally 
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A^ 3 )( g , s ) = e^ 3, ^ <3) l piD ( s ) 



exp 



P3 

1 



(22) 



dr Ap^ 3 ' 3) (r,0;r + s) 



1 - exp / ds' AV(r,s' 



-^y|drAp( 3 ' 4) (r,0;r + s) 1 - exp (-^ J' ds'AV(r, s')) 



Equations (|14|), (19), and (]2^) are the final results of the present theory for the FSE in 3 He- 4 He mixtures. They 
constitute the generalization of the Gersch-Rodriguez formalism to a mixture with special emphasis in the difficulties 
arising from Fermi statistics. Apart from the interatomic potential, very well-known in helium, the microscopic inputs 
that are required are the one- and two-body density matrices, both in the boson-boson and the fermion-boson cases. 

To conclude this section, we define the Compton profiles of each component in the mixture. Contrarily to what 
happens in a pure phase, the total response of the mixture can not be written in terms of a single scaling variable Y. 
Each individual profile is naturally given in its own scaling variable Y a = m a uj/q — q/2. Thus, 



1 r 

jia) ^Y a ) = ^j d S e-^ s S^(q,s) , 
which after introducing the explicit expressions for S( a \q, s) becomes 

(q, Y a )= dY a J( Q ) (Y a )R^ (q, Y a ) + AJ (3) (q, Y a )S a3 . 



(23) 



(24) 



In this equation, AJ^ 3 ) derives from AS^ 3 ) and the IA responses j( a '(Y a ) are directly related to the momentum 
distributions n^ a ' (k) 



J {a \Y a ) =n S(Y i )6 ai + 



4ir 2 p a 



dppn^ a \p) , 



(25) 



no being the 4 He condensate fraction, and v a the spin degeneracy of each component (2/3 = 2, v± = 1). Notice that 
the first term in Eq. (|24| ) contains the explicit contribution uqR^ (q, Y a )S a 4 arising from the condensate. 



III. ENERGY WEIGHTED SUM RULES AT HIGH MOMENTUM TRANSFER 



Energy-weighted sum rules provide an useful tool to analyze the properties of S(q,u>). In spite of the fact that 
the knowledge of a small set of energy moments usually is not enough to completc]jz_x;haracterize the response, the 
method has proved its usefulness in the analysis of scattering on quantum fluids SELa Moreover, from a theoretical 
viewpoint the comparison between the sum rules derived from an approximate theory and the exact ones shed light 
on the accuracy of that approach. In the high-g limit, the response is fully incoherent and therefore we discuss only 
the incoherent sum rules 



/oo 1 in 

Ju^S^u) = __<4 Q C V) [t=0 



Considering 



S^(q,t) = (e-^<"'e mt e lcl <"e-* Ht ) , 



(26) 



(27) 



and applying to the three rightmost operators in (|27|) the Baker-Campbell-Hausdorff formula one arrives to the 
following expansion in terms of it: 



+ ^t) 3 (e-^[H, [H, [H t e^ a) ] 



(28) 
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From Eqs. ( |2q ) and ([28|), one easily identifies the lowest-order sum rules: 

"4?iL(9) = 1 ( 29 ) 

™[% c (q) = W, ^"'l) = ^ (3°) 

mii c (q) = (e-^lH, [H,e^ a) ])) = ( £-) + *^t Q (31) 



2m a J 3 2m a 



,(«) 

l 3,inc 



-C«) r „ r „ ,„ ;„.„(«) 



(q) = (e-^ [H,[H,[H,e^ }}}) = (32) 



2m r 



+ 4(^-)% Q + ^p/ d r 5 (-) W(q .V)V W 



All four moments can be readily calculated from the interatomic pair potential V(r), the kinetic energies per particle i Q , 

and the two-body radial distribution function between pairs of atoms of the same kind g( a ' a \r). [ nc (q) is identical 

to the total m < f\q), also known as the f-sum rule, whereas the other three coincide with the leading contribution to 

the total sum rules m„ (q) at high q. 

In the limit q — > oo the IA is expected to be the dominant term. This feature may be analyzed using the sum-rules 
methodology. Starting from the IA response 

SiaW) = 7^ frfW«)(^(M!-f , (33) 
{2nYp a J \ 2m a 2m a J 

one can calculate the first energy moments from basic properties of the momentum distributions. The results are: 

m$ A (q) = 1 (34) 

(36) 

41W=f^-V + 4f^V*a (37) 



When the IA sum rules are compared with the incoherent results (^9 J30 3^^), one realizes that the first three 
moments are exhausted by IA. The leading order terms in q in the 777.3 sum rule are also reproduced by the IA but 
the term with g( a > a )(r) is not recovered. 

The variable that naturally emerges in the 1 jq expansion of the response of the mixture is the West scaling variable 
Y a . It is therefore also useful to consider the Ya-weighted sum rules of jM(q,Y a ) 



(q) = / dY a Y^J^ (q, Y a ) . (38) 

J —OO 

The first Y a incoherent sum rules are 

(39) 
(40) 

(41) 

y dr</ Q ' Q )(r)(q- V) 2 V{r) . (42) 
In the IA, M (a) , M[ a \ and M 2 (a) coincide with the incoherent sum rules (|9|j4^Q but M^f A = 0. The latter result 





(?) = 


1 


M[ a 


(«) = 





Mt 


(?) = 


2m a 
3 


Mt 


(«) = 


777 a p 

2q 



is in fact general for all the odd Fa-weighted sum rules in the IA due to the symmetry of the IA response around 

r 4 = o. 



G 



In a FSE convolutive theory, as the Gersch- Rodriguez one, it is easy to extract the first sum rules of R(q, Y). From 
the total and the IA sum rules, the use of the algebraic relation 



k 



M k (q) = ( k ) M itIA (q) Mk_ i)R (g) (43) 



4=0 



allows for the extraction of Mj^(q): 

M , R (q) = 1 (44) 
M hK {q) = (45) 
M 2tR (q) = (46) 

M 3 , R (q) = ^p J drg(r)(q-V) 2 V(r) . (47) 

It can be proved that in the Gersch-Rodriguez prescription, the four moments ( ^4||45|j46l , ^7| ) are exactly fulfilled.0 It 
is worth noticing that Ma R (q) is satisfied if and only if a realistic two-body density matrix is used in the calculation 
oUl(q, Y). 

The theory proposed for 3 He in the mixture (Sect. II) predicts a response which is a sum of a convolution product 
plus a correction term AS^ 3 \ The function R^(q,Y 3 ) satisfies M ,r(<7), Mi tR (q), and M 2j r((?) but not M 3 ^ R (q) 
because the convolutive term relies on a boson-boson approximation. Concerning the additive term AS^ 3 -*, it is 
straightforward to verify that their three first moments are strictly zero whereas M^(q) contains corrections to the 
boson-boson g( 3 ' a \r) functions assumed in M 3 _ R (q). 



IV. RESULTS 



The generalization of the Gersch-Rodriguez formalism to the 3 He- 4 He mixture presented in Sect. II requires from 
the knowledge of microscopic ground-state properties of the system. In the present work, the necessary input has 
been obtained using the FHNC/HNC theory.OLJ The variational wave function is written as 

* = F$ 0) (48) 

with F an operator that incorporates the dynamical correlations induced by the interatomic potential, and $0 a model 
wave function that introduces the right quantum statistics of each component. $0 is considered a constant for bosons 
and a Slater determinant for fermions. In the Jastrow approximation, the correlation factor F is given by 

F = F ^UI[f2 a ' m (r l3 ). (49) 

a</3 i<j 

A significant improvement in. the, variational description of helium is achieved when three-body correlations are 
included in the wave functionBj'EZI In this case, 

f = f, t = nn/ 2 M (ni) n n ft^^^r^ . ^ 

a<0i<j a<f3<~fi<j<k 

The isotopic character of the mixture makes the interatomic potential between the different pairs of particles be 
the same. Therefore, the correlation factors f^'^ and f^^^ can be considered to first order as independent of the 
indexes a, (3, 7. That approach, known as average-correlation approximation (ACA),E3 has been assumed throughout 
this work. DMC calculations of 3 He- 4 He mixtureald have estimated that the influence of the ACA in the momentum 
distributions is less than 5 %. 

The dynamic structure function of the mixture has been studied at 3 He concentrations x = 0.066 and x = 0.095 
that, following the experimental isobara P = 0, correspond to the total densities p = 0.3582 c~ 3 and p = 0.3554 cr~ 3 
(a = 2.556 A), respectively. Notice the decrease of p when x increases; in pure 4 He, p = 0.3648 a~ 3 . In Table I, 
results for the 4 He condensate fraction and kinetic energies per particle are reported in J and JT approximations. 
The condensate fraction increases with x whereas the kinetic energies t a decrease, both effects mainly due to the 
diminution of the density. Results for pure 4 He in the JT approximation (the one used hereafter) compare favourably 
with DMC data from Ref. p^4po = 0.084, £4 = 14.3 K), and the decrease of no with x is in agreement with the change 
in no estimated using DMGE2I 
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A. Impulse Approximation 



One of the characteristic properties of the IA in a pure system is its ^-scaling. In this approximation, the response 
is usually written as the Compton profile J(Y). However, global scaling is lost in the mixture due to the different 
mass of the two helium isotopes. The individual Compton profiles j( a \Y a ) must be written in terms of its own Y a 
variable. 

Results for J {a) (Y a ) at x = 0.095 are shown in Fig. 1. The different statistics of 4 He and 3 He are clearly visualized 
in their respective momentum distributions, and therefore also in the Compton profiles. In (Y±), a delta singularity 
of strength no located at I4 = (not shown in the figure) emerges on top of the background, whereas in J^ 3 \Y 3 ) 
the Fermi statistics is reflected in the kinks at Y 3 = ±kp produced by the gap of n^(k) at k = kp. The large 
\Y a \ behavior of both responses is more similar and is entirely dominated by the tails of the respective momentum 
distributions. 

The dwnamic structure function of the mixture suggests the definition of a total generalized Compton profiles 
J(q,Y a )B In the IA, 

J(q,Y a ) = — 1 — Si A (q,uj) , (51) 

cr a (d a3 x + d a4 {l —x)) m a 



with 



5iA(g,w) = <74(1 - x) S$(q, uj) + <r 3 x S$(q,uj) . (52) 



Notice that the definition (|T]) is different for each Y a . In order not to overload the notation, the introduction of a 
new labelling in J(q, Y a ) has been omitted. In terms of Y±, and introducing the single Compton profiles J^{Y a ), 

J{q,Y±) = j( 4 >(Y 4 ) + — ^ ^ J^{Y 3 (Y,)) , (53) 
cr 4 (l — x) rri4 



with 



Y 3 (Y 4 ) = ^Y 4 -Ul-^) ■ (54) 



772,4 2 \ 7774 

Equivalently, one can express the total generalized Compton profile as a function of Y 3 



J {q ,Y 3 ) = ^^l^±J^(Y 4 (Y 3 )) + J(3)(r 3 ) , (55) 
cr 3 x m 3 



with 



W^S-^P-l ■ (56) 
m 3 2 \m 3 ) 

The choice of the scaling variable Y a undoubtedly determines some trends of the response. If I4 is used, the 4 He 
peak is centered at Y4 — and the 3 He peak shifts to I4 = (7714/7713 — I) q/2 ~ q/6. On the other side, if Y 3 is the 
choice the 3 He peak is centered at Y 3 = and the 4 He one moves to Y 3 — (7713/7714 — 1) q/2 ~ —q/8. In addition, and 
disregarding cross sections and concentration factors, the 3 He peak is reduced by a factor ra 3 jra 4 when the response 
is expressed in terms of Y4. By the same token, the 4 He peak is enhanced by a factor 7714/7713 when the response is 
written as a function of Y 3 . 

In Fig. 2, the IA responses for the mixture at two different 3 He concentrations are shown. They correspond to 
a momentum transfer q = 23.1 A^ 1 and have been obtained from 7i^ Q )(fc) calculated at the JT approximation level. 
The differences between both curves are due to the concentration factors rather than to the differences between the 
momentum distributions involved. 



B. Final State Effects 



The theory of FSE in 3 He- 4 He mixtures developed in Sect. II requires from the knowledge of the three correcting 
functions R {4) {q, s) @, i? (3) (q, s) @, and AS" (3) (q, s) @ (s = tq/m a ). These three functions are complex 
with real and imaginary parts that are, respectively, even and odd functions under the change s — * — s. The latter 
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is a consequence of the symmetry properties of the two-body density matrices and of the central character of the 
interatomic potential. The Fourier transforms of the real and imaginary parts generate, respectively, the even and 
odd components of R^ a \q, Y a ) and AS^ 3 )(q, Y 3 ), which are all real. 

In Fig. 3, the real and imaginary parts of R( a \q,s) corresponding to a x = 0.095 mixture are shown. In spite 
of the fact that R^(q,s) is calculated for the real mixture and R^ 3 \q,s) for the boson-boson one, the differences 
between the two functions are rather small. Actually, those differences are mainly attributable to the low 3 He density 
in the mixture that makes the contributions of the Fermi statistics very small. In fact, the differences shown in Fig. 3 
between R^ A \q, s) and R^ 3 \q, s) are essentially due to the different mass of the two isotopes, which factorizes in the 



integral of the interatomic potentials (see Eqs.l4[l9|). 

The real and imaginary parts of the additive term AS^(q, s) are shown in Fig. 4 at the two 3 He concentrations 
studied. The behavior of AS^ 3 \q, s) is remarkably different from the behavior of the FSE broadening functions 
R( a \q, s), presenting oscillating tails that slowly fall to zero with increasing x. The function AS^ 3 '(q, s) incorporates 
on the 3 He response all the Fermi corrections which are not contained in R^ 3 '(q, s). In a dilute Fermi liquid, as 3 He in 
the mixture, those contributions are characterized by the behavior of l(kpr) and l 2 (kpr), l(z) — 3/ z 3 (sinz — z cos z) 
being the Slater function. 

i?( 4 )(<7, Y4) and R^ 3 \q, Y3) are compared in Fig. 5 at x — 0.095 and q = 23.1 A -1 . The shape of both functions 
looks very much the same: a dominant central peak and small oscillating tails that vanish with \Y a \. The figure 
also shows that at a given concentration the central peak of R^(q, I3) is slightly higher and narrower than the one 
of i?( 4 )(<?,r 4 ), an effect once again due to the different mass of the two isotopes. Therefore, at a fixed momentum 
transfer q FSE in 4 He are expected to be smaller than in 4 He. In the scale used in Fig. 5, the R^(q, Y a ) functions 
at x = 0.066 would be hardly distinguishable from the ones at x = 0.095. 

The Compton profile A (q, Y 3 ), derived from the Fourier transform of AS^ 3 ) (q, s), is shown in Fig. 6 at the two 
x values considered. A (q,Y 3 ) presents a central peak and two minima close to Y3 = ±fci?- The absolute value 
of this function is small compared to both R^(q : Y 3 ) and the IA response J^(Y 3 ) ( p5| ) but manifests a sizeable 
dependence on the 3 He concentration. This feature is patent in Fig. 6, where one can see how the contribution 
of AJ^(q,Y 3 ) increases with x. This is an expected result taking into account that in the current approximation 
Aj( 3 )(q, Yi) incorporates all the Fermi effects to the 3 Hc FSE function. 

According to the theory developed in Sect. II, the 4 He response in the mixture, J^'(q, Y4) is the sum of two terms: 
the non-condensate part of the IA convoluted with R( 4 \q,Y4), and n^R^ (q, I4), which is the contribution of the 
condensate once broadened by FSE. The different terms contributing to the final response are separately shown in 
Fig. 7. The correction driven by no is by far the largest one. In spite of the small value of no, the broadening of 
the condensate term, which transforms the delta singularity predicted by the IA into a function of finite height and 
width, unambiguously produces non-negligible FSE in the 4 He peak. 

The obvious lack of a condensate fraction in the 3 He component reduces the quantitative relevance of its FSE. 
The 3 He FSE correcting functions and the corresponding IA response, are compared in Fig. 8 at x — 0.095. The 
convolution of the IA with R^ 3 \q,Y 3 ) produces a slight quenching of J^ 3 \q,Y 3 ) around the peak and a complete 
smoothing of the discontinuity in the derivative of J^(Y}) at Y 3 = ±kp. The contribution of A (q,Y 3 ) is rather 
small but restores to some extent the change in the derivative around kp . 



C. Theory vs. Experiment 

Scattering experiments suffer from instrumental resolution effects (IRE) that tend to smooth the detailed structure 
of the dynamic structure function. Any comparison between theory and experiment have therefore to include in the 
analysis the IRE contributions. From the theoretical side, it would be desirable to remove the IRE from the data 
to allow for a direct comparison. This ^process would imply a deconvolution procedure that is known to be highly 
unstable. As suggested by Sokol et a/.,Ea it is better to convolute the theoretical prediction with the IRE function. 
I [a) {q,Y a ), and then to compare the result with the experimental data. The functions I^ a \q, Y a ) provided by Sokolc2l 
are reported in Fig. 9. As one can see, at q = 23.1 A -1 the IRE corrections are of the same order of the FSE functions 
R( a '(q,Y a ), and in fact their magnitude significantly increases with q. The IRE functions for the mixture (Fig. 9) 
present a small shift of their maximum to negative Y values, a feature that makes the peak of the total response be 
slightly moved in the same direction. 

In Fig. 10, the generalized Compton profile J(q : Yi) (including both the IRE and FSE) is compared with the 
scattering data of Wang and SokoLu Those measurements were carried out in a x = 0.095 mixture at T = 1.4 K 
and a momentum transfer q = 23.1 A -1 . The analysis of the experimental data led the authors to estimate the 4 He 
condensate fraction and the single-particle kinetic energies of both species. In Ref. ||, a value no = 0.18 ± 0.03, and 
kinetic energies £4 = 13 ±3 K and t 3 = 11 ±3 K are reported. That work, and an independent measurement performed 



9 



by Azuah et al.f^i agree in the values of the kinetic energies and in their dependence with the 3 He concentration. Both, 
analysis coincide in a decrease in t± with x and a more surprising constancy of t 3 along x. Microscopic ealculationsEj 
of those quantities only agree with the experimental result of t±(ft). Several independent calculations,™^ including 
the present one, suggest smaller values of no (no ~ 0.10) and larger values of £3 (£3 ~ 18 K), in clear disagreement 
with the experimental estimations. 

Let us turn to Fig. 10 with the comparison between the theoretical and experimental responses. The theoretical 



result, constructed using Eqs. (HH), but replacing the IA J {a \Y a ) with the final responses jW(<?,y a ), shows 
sizeable differences with respect to the experimental data and a lack of strength below the two peaks. In order 
to clarify the origin of such a large discrepancy, we have compared the Mo and Mi sum rules obtained by direct 
integration of the experimental J(q, I4) with the theoretical results (Sect. III). That check has shown that the Mo 
and Mi values obtained from the two procedures are not compatible. Our conclusion is that the reported experimental 
Compton profiles are probably written in in a different way than in Eq. (fill). In fact, after the analysis of different 
possibilities, we have verified that if one defines the response in the form 

J(q,Y A ) = J (4 \q,Y 4 ) + °* X J^(q,Y 3 (Y 4 )) (57) 
0-4(1 — x) 

or 

J(q,Y 3 ) = ~ X) jW( g ,y 4 (y 3 )) + jW( q> Y 3 ) , (58) 

<T 3 X 

the agreement in both sum rules is recovered. By moving our results to those modified Compton profiles J(q,Y a ), 
the agreement between theory and experiment improves significantly but only to what concerns the 3 He peak. Notice 
that the 4 He peak is not modified when going from J(q, Y4) to J(q, ¥4), and that a significant difference in the height 
of the peak still remains. 

The missing strength of the theoretical 4 He peak with respect to the experimental data could justify the difference 
between the theoretical and experimental values of no. However, the .present variational momentum distribution 
predicts no values that are indistinguishable from a DMC estimation.llj Therefore, this difference should not be 
attributed to inaccuracies of our n^ a '(k) but rather to an intriguing gap between theory and experiment. At this 
point, it is worth considering the difficulties the experimentalists have to face to extract no and t a from the measured 
data. On the one hand, experience in the pure 4 He response has shown that different momentum distributions (with 
different no's) can be accurately fitted to the data. On the other, the kinetic energy per particle is derived from the Y£ 
sum rule whose estimation is highly influenced by the tails of the response. Those tails cannot be accurately resolved 
due to the noise of the data, and thus the prediction of t a appears relatively uncertain. That is even more pronounced 
in the 3 He peak because the strong interaction with 4 He causes n^(k) to present non-negligible occupations up to 
large k values. 

The influence of no and t a on the momentum distribution, and hence on the response, can be roughly estimated 
from the behavior of the one-body density matrix. In a simple approximation, one can perform a cumulant expansion 
of Pi (r) and relate the lowest order cumulants to the lowest order sum rules of n^ a \k). Introducing an expansion 
parameter A, 



-^ 4) (Ar)-n = e ^ +A2 ^ + - = (l-no)-A 2 ((r. Pl ) 2 ) + ..- . (59) 
Pi 



Taking into account that 



and considering A = 1, 



2\ Im^r 2 



((r-pi) 2 ) = ^—U 7 (60) 



— p^ 4) (r) = n + (1 - no) exp 

P4 



2 m^r 



2 



3 (1 - n ) 



u + 



(61) 



Equation (|6|) can then be used to relate pf\r) to a new one-body density matrix Pi (r) with slightly different 
values no and I4 



— Pi 4) ( r ) = n + ( — p{ 4) (r) -n ) ( \ — — ) exp 
P4 \PA J \l-no ' 



Im^r 2 ( ?4 ti 



1 - n 1 - n 



(62) 
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In this way, the perturbed Pi (r) and n*- 4 -* (fc) preserve their normalization and allows one to go beyond a simple uq 
re-scaling. Using this method, we have studied the effect of changing no and t± on the 4 He response. In Fig. 11, the 
results corresponding to i) n = 0.14, t 4 = 13.9 K, and it) n = 0.10, t 4 = 13.0 K are shown. As one can see, both 
slight changes in the theoretical response lead to a nice agreement with the experimental data. Consequently, such a 
large value of no (ng Xpt = 0.18) does not seem to be required in order to reproduce the additional strength observed 
below the 4 He peak. The re-scaling (|6^) shows that a small decrease in the kinetic energy enhances the central peak 
in the same form an increase of the condensate fraction docs. 



V. SUMMARY AND CONCLUSIONS 



A generalized Gersch- Rodriguez formalism has been applied to study the dynamic structure function of the 3 He- 4 He 
mixture at high momentum transfer. The Fermi character of 3 He forbids a straightforward generalization of most 
FSE theories used in bosonic systems, a problem that has been overcome in an approximate way. The approximations 
assumed are however expected to include the leading Fermi contributions to the FSE, at least in the mixture where 
the 3 He partial density is very small. 

The theoretical response obtained shows significant differences with scattering data in both the 4 He and the 3 He 
peaks. However, a sum-rules analysis of the experimental response has shown some inconsistencies. Redefining the 
total response, it is possible to reach agreement between the theoretical and the numerical values of the first-order 
sum rules. If the theoretical response is changed in the same way, the agreement is much better. Nevertheless, the 
4 He peak is not modified by this redefinition (written as a function of Y4) and an intriguing sizeable difference in 
its strength subsists. From the theoretical side, several arguments may be argued trying to explain the observed 
discrepancies. The first uncertainty could be attributed to the use of a Gersch-Rodriguez theory to account for the 
FSE. In our opinion, that criticism has probably no sense because we have verified that, at similar momentum transfer, 
the experimental response of pure 4 He is fully recovered with the GR theoryO Assuming therefore that the theoretical 
framework is able to describe the high-q response of the mixture, one could be led to argue that the approximate 
microscopic inputs of the theory are not accurate enough. That argument was put forward in Ref. |^ to explain the 
differences in hq and £3. One of the main criticisms was the use of the AC A, which they claimed could be tocujestrictive 
to allow for a reduction of t§ towards a value closer to the experimental one. However, a DMC calculation^!] in which 
the ACA is not present, has proved that only a diminution of ~ 0.5 K in t 3 is obtained. Concerning the condensate 
fraction value, our variational theory predicts a slight increase of no with x. This increase, which is mainly due to the 
decrease of the equilibrium density when x grows, is nevertheless much smaller then the one that would be required 
to reproduce the experimental prediction. Our results for no & r e again in an overall agreement with the nearly exact 
DMC calculation of Ref. [l0|. 

In summary, we would like to emphasize that there exists theoretical agreement on the values of no and £3 for 
mixtures, but these values are quite far from the experimental estimations. Additional scattering measurements on 
the 3 He- 4 He mixture are necessary to solve the puzzle. 
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TABLE I. Condensate fraction and kinetic energies as a function of x. At each 3 He concentration x the first row corresponds 
to the J approximation and the second one to the JT one. 
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FIG. 1. Compton profiles of 4 He (left) and 3 He (right), both in JT (solid line) and J (dashed line) approximations for 
2=0.095. 
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FIG. 2. Generalized Compton profiles in IA at x = 0.066 (left) and x — 0.095 (right). 
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FIG. 3. Real and Imaginary parts of R^fas) (solid line) and R^(q, s) (dashed line) at q = 23.1 A 1 for the x = 0.095 
mixture. 
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FIG. 4. AS (3) (g, s) at g = 23.1 A -1 and for mixtures at x = 0.095 and x = 0.066 (solid and dashed lines). 



0.0060 



0.0040 



& 0.0020 



CO 

< 

p§ 0.0000 



-0.0020 




0.0010 



0.0005 g 
> 

0.0000 ^ 

CO 

I 

S 

-0.0005 3 



20 40 
s(hA) 



20 40 
s(hA) 



60 



-0.0010 



14 



FIG. 5. Comparison between (q, Y4) and R^ 3 \q, Y3) at q = 23.1 A 1 and for x = 0.095 (solid and dashed lines, respec- 
tively). Notice that different Y a variables are used to depict each function. 
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FIG. 6. The 3 He additive correcting term at q = 23.1 A 1 for x = 0.095 and x = 0.066 mixtures (solid and dashed lines, 
respectively) . 
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FIG. 7. The different contributions to the 4 He response at x = 0.095 and q = 23.1 A 1 . Dotted line: 4 He Compton profile, 
dashed line: the same convoluted with R( 4 \q, Y 4 ), long-dashed line: n R^ (q, Y 4 ), solid line: total 4 He response. 




FIG. 9. 4 He (left) and 3 He (right) instrumental resolution functions at x = 0.095 and q = 23.1 A 




FIG. 11. The x = 0.095 experimental data of Wang and Sokotl compared to the response obtained from an alternative 
~p^\r) with no = 0.14 and I4 — £4 = 13.9 K (solid line) -left panel-, and with no = 0.10 and I4 — 13.0 K (solid line) -right 
panel-. 
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